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1 Introduction 

New CO observations at the center of cooling flow in Perseus cluster [4] show 
a clear correlation of the molecular gas with the previously detected H-cv fil- 
aments [1]. In this poster, we present high resolution multi-phase simulations 
of the Perseus Cluster, taking into account the AGN feedback in form of 
hot buoyant bubbles. These simulations show that significant amount of gas 
can cool far from the center. The AGN feedback provides some heating, but 
also trigger the hot gas compression, that favours cooling (positive feedback) , 
even at high radius (R > 30 kpc) . The cooled gas flows into the cluster core 
forming the observed filaments. 

2 Initial Conditions and gas physics 

Our cluster model is designed to fit the Perseus X-ray data [5]. The total 
mass distribution profile follows a pseudo-isothermal sphere with a central 
electronic density of 5 x 10~ 2 and a core radius of 40 kpc. The mass dis- 
tribution is truncated at 2Mpc and the total mass is 5.5 x 1O 14 M0. The 
gas corresponding to 15% of the total mass has an initial temperature of 
2.8 x 10 7 K. 

The bulk of the intra-cluster gas has at high temperature (10 7 K) is well 
modeled by an ideal gas with adiabatic index of 5/3. This hot phase has 
been computed using and SPH technique. For lower temperatures, SPH is 
not suited to render the dumpiness of the gas. Below 10 4 K, the gas is treated 
as semi-collisional, using a sticky particle technique [2]. Star formation may 
occur in this cold phase, following a Schmidt law. The cooling of the hot gas 
is computed using the standard cooling function [8] , assuming an abundance 
of one third of the solar one. 

Following [6], the AGN feedback is modeled by injecting energy in hot 
intra-cluster gas bubbles which are then driven buoyantly at higher radius. 
Bubbles are defined by their position, diameter, temperature, over pressure 
and angle. They are generated by symmetric pairs, in average each 200 Myr. 

The simulations have been run using a modified version of the Gadget2 
Tree-SPH code (Springel, 2005), including cooling, star formation and multi- 
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phase SPH/Sticky gas. In each simulations the gas is represented by 4'194'304 
SPH particles. The dark matter is modeled by a fixed outer potential. 

3 Results 

Fig. 1. shows the evolution of an isolated bubble. The bubble has an initial 
temperature of 10 8 K. It reaches 30 kpc and takes the characteristic mush- 
room shape as observed in H-a and X-ray [3] . A substantial amount of cooler 
gas (about 1/2 of the bubble mass) is dragged by the rising bubble and forms 
the trunk of the mushroom. Between t = 300 — 500 Myr, this gas falls back 
into the cluster center, while the head of the mushroom stops and dissipate 
into the intra-cluster medium. 
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Fig. 1. Temperature map of an isolated bubble. From 10 4 K to 2 x 10 8 K. The box 
dimension is 200 x 200 kpc. 



3.1 Global evolution and cold gas formation 
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Fig. 2. Global evolution of the cluster temperature during 4Gyr. The box dimen- 
sion is 400 x 400 kpc. 



In addition to slowing down the cooling flow at the center of the cluster, 
our simulation also show how AGN feedback may trigger cold gas formation 
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at high radius (R > 30 kpc) . Physical conditions for efficiently cool intra- 
cluster gas occurs either when the gas is less than 10 7 K or when its density 
is sufficiently high. Buoyant bubbles are responsible of strong inhomogeneities 
in temperature (see Fig. 2) as well as in density of the intra-cluster medium. 
In Fig. 3 left, we show how cooler falling gas (trunk of a bubble, see Fig. 2) is 
trapped between an old and a new rising bubble and is compressed to reach a 
state where its cooling time is sufficiently short to let the gas becomes cool in 
a fraction of Gyr. In Fig. 3 right, the cluster is seen 300 Myr later. The short 
cooling time gas of Fig. 3 left, has now cooled down below 10 4 K. As it is not 
supported by the hot gas pressure, it falls radially to the center, forming a 
filament like structure (50 to 100 kpc) of mass 1.5 x 10 8 Mq. In the filament, 
the gas density is not high enough to form stars. 




Fig. 3. Left : Gas temperature at t — 2720 Myr. The white contours pointed by 
black arrows indicate the position of the gas having a short cooling time r c < l.l Gyr 
that will be transformed into cold gas. The upper left box corresponds to the region 
zoomed in the right part of the figure. Right : Gas temperature between t — 2950 
and t = 3100 Myr, from 10 4 to 2 x 10 8 K . The small black dots pointed by the 
black arrow represents cold gas (T < 10 4 K) falling into the cluster center. 
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